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pH ( Aims. We study the synchrotron flaring behaviour of the blazar 3C 279 based on an extensive dataset covering 10 years of monitoring at 19 
Q | different frequencies in the radio-to-optical range. 

Methods. The properties of a typical outburst are derived from the observations by decomposing the 19 lightcurves into a series of self-similar 
events. This analysis is achieved by fitting all data simultaneously to a succession of outbursts defined according to the shock-in-jet model of 
Marscher& Gear (1985). 

Results. We compare the derived properties of the synchrotron outbursts in 3C 279 to those obtained with a similar method for the quasar 
3C 273. It is argued that differences in the flaring behaviour of these two sources are intrinsic to the sources themselves rather than being 
due to orientation effects. We also compare the start times and flux densities of our modelled outbursts with those measured from radio 
k> components identified in Very Long Baseline Interferometry (VLBI) images. We find VLBI counterparts for most of our model outbursts, 
*j_J . although some high-frequency peaking events are not seen in the radio maps. Finally, we study the link between the appearance of a new 
C3 ' synchrotron component and the EGRET gamma-ray state of the source at 10 different epochs. We find that an early-stage shock component 
is always present during high gamma-ray states, while in low gamma-ray states the time since the onset of the last synchrotron outburst is 
significantly longer. This statistically significant correlation supports the idea that gamma-ray flares are associated with the early stages of shock 
components propagating in the jet. We note, however, that the shock wave is already beyond the broad line region during the gamma-ray flaring. 
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1. Introduction 

It is generally accepted that radio outbursts in active galactic 
nuclei (AGN) are triggered by growing shocks in a relativis- 
tic jet. There have been many efforts to identify the individ- 
ual components in radio to submillimeter (submm) lightcurves. 
Litchfield et al. (1995) and Stevens et al. (1995, 1996, 1998) 
managed to follow the early evolution of individual outbursts 
by subtracting the quiescent contribution (assumed to be con- 
stant) from the total spectrum. Valtaoja et al. (1999) identified 
outbursts by decomposing the variations in millimeter (mm) 
and centimeter lightcurves into exponential flares. Tiirler et al. 
(1999) measured the properties of synchrotron outbursts by de- 
Send offprint requests to: EJ. Lindfors email:elilin@utu.fi 



composing the radio-to-submillimeter wavelength lightcurves 
into a series of self-similar flaring events. The method was 
further developed by Tiirler et al. (2000) who used the shock 
model of Marscher & Gear (1985) to describe both the aver- 
age evolution of the outbursts and their individual characteris- 
tics. This model was found to provide a good description of 
the lightcurves of 3C 273. In this work, we decompose the 
lightcurves of a second source, 3C 279. Since our adopted 
methodology is analogous to that used by Tiirler et al. (2000) 
the results are directly comparable to those found for 3C 273. 

While in previous studies the outbursts were identified in 
the radio-to-submillimeter regime, we here follow the syn- 
chrotron spectrum up to infrared and optical frequencies. 
Searches for correlations in the radio-to-optical emission from 
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AGN have been conducted since the 1970s. One of the most re- 
cent works was presented in Hanski et al. (2002). In their sam- 
ple of 20 AGNs they found a clear radio-to-optical correlation 
in seven cases, a possible correlation in six cases and no corre- 
lation for the remainder. This agrees well with the findings of 
previous studies: there appears to be some kind of correlation 
between the radio and optical emission but this correlation is 
not seen in all sources at all epochs. The general trend seems to 
be that all radio outbursts are accompanied by optical outbursts 
but conversely not all optical outbursts have radio counterparts. 

The source 3C 279 is one of the brightest and most variable 
blazars at radio wavelengths. As such, it is monitored regularly 
in all radio wavebands from the centimeter to the submillime- 
ter domain. In the optical, the historical lightcurve shows varia- 
tions with a typical amplitude of about 2 magnitudes, but reach- 
ing 8 magnitudes during flares (Webb et al. 1990). The weak 
blue bump of 3C 279 allows us to follow the synchrotron spec- 
trum up to infrared and optical wavelengths. Indeed, 3C 279 is 
one of the few objects for which a clear correlation is found be- 
tween the radio and optical wavebands (Tornikoski et al. 1994). 

In this paper we present results of a multifrequency 
lightcurve decomposition of 3C 279. The method used is based 
on Tiirler et al. (1999; 2000) with some modifications for the 
specific case of 3C 279, including an extension of the analy- 
sis to infrared and optical wavelengths. The average evolution 
of outbursts is compared to that of 3C 273. We also compare 
our results with (1) VLBI maps from Wehrle et al. (2001) and 
Jorstad et al. (2004) to study the connection between the model 
outbursts and VLBI components, and (2) gamma-ray data ob- 
tained by Hartman et al. (2001) to study the connection be- 
tween the synchrotron spectra and the gamma-ray state. 

2. Data 

To decompose the radio-to-optical lightcurves of 3C 279 we 
used data from five radio and seven mm/submm wavelengths 
as well as four infrared and two optical wavebands. The data 
sample extends from 1989.0 to 1999.5. 

The 4.8, 8.0 and 14.5 GHz data are from the University 
of Michigan Radio Astronomy Observatory (UMRAO). Some 
of these data have not previously been published. The 22 
and 37 GHz data and some of the 90 GHz data are from the 
Metsahovi Radio Observatory. These data were published by 
Terasranta et al. (1998, 2004). We have also used 90 and 
230 GHz (1.3 mm) data from the Swedish-ESO Sub-millimeter 
Telescope (SEST). SEST data taken prior to 1994.5 were pub- 
lished by Tornikoski et al. (1996) while more recent data are 
published here for the first time. From the Institut de Radio 
Astronomie Millimetrique (IRAM) we have data at 90, 150 
and 230 GHz (Steppe et al. 1993; Reuter et al. 1997). We also 
used data at 2.0, 1.3, 1.1, 0.85, 0.8, 0.45 and 0.35 mm from the 
James Clerk Maxwell Telescope (JCMT). The 0.85 mm dataset 
was published by Robson et al. (2001) while data in the other 
wavebands are either from Stevens et al. (1994) or are previ- 
ously unpublished (post 1993). The dataset of Stevens et al. 
(1994) also includes the four infrared wavebands used in this 
work. 



The optical R- and V-band lightcurves are collected from 
the literature, but also contain previously unpublished data. We 
have used the data from Maraschi et al. (1994), Hartman et 
al. (1996), Villata et al. (1997), Katajainen et al. (2000) and 
Hartman et al. (2001). The previously unpublished data are 
from the Kungliga Vetenskapsakademien (KVA) telescope and 
the Nordic Optical Telescope (NOT) as well as from the obser- 
vatories of Perugia and Torino. 

3. Model and Method 

In this paper we study the multifrequency lightcurves of 3C 279 
in the context of the shock model of Marscher & Gear (1985). 
As the shock propagates downstream in the relativistic jet it 
evolves through three stages. In the initial (growth) stage, when 
inverse Compton losses predominate, the synchrotron self- 
absorption turnover frequency decreases and the turnover flux 
density increases. In the second (plateau) stage, synchrotron 
losses dominate and the turnover frequency decreases while 
the turnover flux density remains roughly constant. During 
the third (decay) stage when adiabatic losses dominate, both 
turnover frequency and turnover flux density decrease. 

Bjornsson & Aslaksen (2000) criticize one of the assump- 
tions of Marscher & Gear (1985) which concerns the rise of the 
peak flux density during the initial Compton stage. Their cor- 
rection of the expression for the energy density of synchrotron 
photons implies a much shallower rise of peak flux density with 
peak frequency at the beginning of the outburst. On the other 
hand, multiple Compton scattering, a process not considered by 
Marscher & Gear (1985), would actually steepen their derived 
relation. There is therefore no compelling theoretical reason to 
use a modified version of the original Marscher & Gear (1985) 
model. 

In Tiirler et al. (2000), a generalisation of this model was 
used to decompose the multifrequency lightcurves of 3C 273. 
We adopt a similar methodology in this work. The shape of 
the emission spectrum behind the shock front is assumed to 
be that of a simple synchrotron spectrum (with electron en- 
ergy distribution N(E) oc KE~ S ) with two spectral breaks, 
namely the low-frequency break, Vh, below which a t hi c k = 5/2 
and the high frequency break, Vb, which steepens a±i n from 
(1 - s)/2 to -s/2. To extend the decomposition into the in- 
frared and optical we introduce an additional high frequency 
exponential cut-off to the synchrotron spectrum corresponding 
to emission from the most energetic electrons. The sharp low- 
and high-frequency spectral breaks of Tiirler et al. (2000) were 
also replaced by more realistic smooth transitions. As in the 
modelling of the micro-quasar GRS 1915+105 (Tiirler et al. 
2004), we introduced a new parameter allowing the ratio Vh/v m 
of the low-frequency spectral break, Vh, and the synchrotron 
self-absorption turnover, v m , to vary with time. 

Another significant change is the addition of an underly- 
ing jet component with constant flux density. This compo- 
nent is assumed to be an inhomogeneous synchrotron source 
with a spectrum defined by a high-frequency break fixed at 
375 GHz - as suggested by the quiescent level studies of 
3C 279 (Litchfield et al. 1995) - and by two free parameters of 
the fit defining the synchrotron self-absorption flux density and 
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Fig. 1. Ten out of nineteen radio-to-optical lightcurves of 3C 279. The points indicate the observed flux density and the solid line 
our best fit which is a sum of the underlying jet (long dashed line) and the fifteen outbursts (dotted lines). The first dotted line 
represents the global decay of all outbursts peaking before 1989.0 



frequency. To be consistent with the lowest infrared measure- 
ments an exponential cut-off was also added to the underlying 
jet spectrum at a somewhat arbitrary frequency of 2-10 5 GHz. 
We note, however, that our conclusions are insensitive to these 
numerical values. 

The specificity of the individual outbursts is modelled by 
varying the scale of their evolution in flux density, frequency 
and time as done in Tiirler et al. (1999). The more physical 
approach adopted by Tiirler et al. (2000) of varying the onset 
values of K (electron energy density normalization), B (mag- 
netic field strength) and D (Doppler factor) was also tried, but 



was found to limit the differences between outbursts too much 
to get an equally good fit. 

The model presented here assumes a constant Doppler fac- 
tor D oc R~ d and a magnetic field B oz R h perpendicular to 
the jet axis as suggested by VLBI polarization measurements 
(Lister et al. 1998, Marscher et al. 2002). The exponents, b and 
d, of the evolution of these quantities with radius, R, of the jet 
are therefore fixed to d = and b = 1 . We tried varying these 
values as well, but this did not result in significantly better fits. 

The resulting model has a total of 77 free parameters to 
fit the 15 outbursts identified in the lightcurves of 3C 279 from 
1989 to 2000. Twelve parameters are used to describe the shape 
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and evolution of the synchrotron spectrum for an average out- 
burst. Three further parameters describe the initial flux density 
decay while two are needed to describe the underlying jet com- 
ponent. The remaining 60 parameters are used to describe the 
start time and the specific characteristics of the 15 different out- 
bursts. The total number of degrees of freedom is 3281. 

4. Results 

The best-fit decomposition of the radio-to-optical lightcurves 
into a series of 15 self-similar outbursts is shown in Fig. 1. It 
has a reduced^- 2 value of 8.72. This relatively high value can 
be considered as acceptable here keeping in mind that the aim 
of the modelling is to derive the average properties of a typical 
outburst and that the specificity of each outburst is only mod- 
elled very crudely in order to minimise the number of free pa- 
rameters. Furthermore, the analytical shock-in-jet model itself 
is of necessity a simplistic idealization of the complex emission 
of actual jets (Gomez 2005) and thus a perfect match of the 
model to the data cannot be expected. At the lowest frequen- 
cies, the fit describes the smooth shape of the lightcurves very 
well. At 90 GHz, the general flux density level of the fit is good, 
but during the biggest outbursts (in 1991, 1993.5, 1994 and 
1996.5) the model flux density remains below the observations. 
The same tendency also applies to all mm-band lightcurves and 
to the optical R- and V-band lightcurves. Part of this fit (until 
1994.0) has been published in Lindfors et al. (2005). 

4.1. Average evolution of the synchrotron spectrum 

The average evolution of an outburst in 3C 279, as defined by 
12 fit parameters, is shown in Fig. 2. In Fig. 2c we see the 
average outburst spectra. It peaks at ~ 150 GHz, which is the 
frequency at which the transition from the Compton stage to 
the synchrotron stage takes place. The peak flux density is 6.6 
Jy. The second stage transition in the average spectra is at ~ 
20 GHz and its flux density is 2.2 Jy. 

The overall shape of the evolution is quite different from the 
one derived for 3C 273 (Tiirler et al. 2000). The most obvious 
difference is the decreasing flux density during the second stage 
of the evolution which is quite different from the typical plateau 
seen in 3C 273. In the model, this difference can be ascribed 
to the high value of the parameter k defining the decrease of 
the normalization K oc R k of the electron energy distribution 
N(E) oc KE S with jet opening radius, R. The derived value of 
k is 4.0 which is far above the value of fe ad = 2{s + 2)/3 = 2.81 
(with s = 2.25) predicted for an adiabatic jet flow. The origi- 
nal Marscher & Gear model adopted only adiabatic compres- 
sion to heat the electrons in the shock. While this seems to 
be a good assumption for 3C 273 (k=3.03), our high A:-value 
found for 3C 279 suggests that the emitting electrons are sub- 
ject to important non-adiabatic cooling processes within the jet. 
This might be due, at least in part, to synchrotron and inverse 
Compton radiative losses in the undisturbed underlying jet. 

As mentioned in Tiirler et al. (2000), the three parameters 
r, k and d are difficult to constrain uniquely by the fit. In par- 
ticular, the effect of k > k a d and of d > are similar, both pro- 
ducing a decreasing flux density during the synchrotron stage 



of the evolution. The high value of k might therefore also indi- 
cate that our assumption of d = is not valid and thus suggest 
that the Doppler factor tends to decrease during the evolution 
of the shock. This could be due to a decelerating jet flow or to 
a geometry in which the jet bends away from the line of sight 
(see Jorstad et al. 2004). 

Another difference is that the flattening of the spectral in- 
dex by Aa,hi„ = +0.5 (predicted to occur at the synchrotron- 
to-adiabatic stage transition) is found to be much less abrupt 
for 3C 279 than for 3C 273, and for 3C 279 it starts much 
earlier during the initial Compton stage. This can be seen 
by the position of the start of the flattening (ff) in Fig. 3, 
which is much earlier than found for 3C 273 (see Fig. 4. of 
Tiirler et al. (2000)). A similar behaviour was also found for 
GRS 1915+105 (Tiirler et al. 2004). 

Apart from k, the values of two other physical parameters 
of the jet are found to be similar to those derived for 3C 273 
(Tiirler et al. 2000). The index, s, of the electron energy dis- 
tribution is found here to be s = 2.25 (s = 2.05 for 3C 273) 
and the opening radius R oc U of the jet with distance, L, along 
the flow is found to be slightly non-linear (r = 0.78) indicating 
that, as for 3C 273 (r = 0.8), the jet opening angle is slowly 
decreasing with distance suggesting that some jet collimation 
process is at work. 

We also note, as was found for GRS 1915+105 (Tiirler et al. 
2004), that there is a tendency for the synchrotron emission of 
the outbursts to start out inhomogeneous but to become homo- 
geneous as the shock evolves with time. This is manifested as a 
decreasing Vh/v m ratio which affects the shape of the spectrum 
at frequencies below the turnover as can be seen in Fig. 2c. 

Another difference is that in 3C 273 the outbursts are very 
short-lived at higher frequencies, while in 3C 279 at all fre- 
quencies above ~ 150 GHz it takes ~ 0.45 years to reach the 
peak flux density. Since the rise time at all such high frequen- 
cies is the same these outbursts peak simultaneously. This is 
shown in Fig. 3. This result, together with the finding that the 
plateau stage is a bit shorter than in 3C 273, might suggest that 
in 3C 279 the magnetic field energy density is relatively lower 
than the electron energy density. 

3C 279 is an archetypical blazar while 3C 273 is not a 
true member of the blazar class. Most conspicuously, it is not 
rapidly variable in the optical band and has low optical polar- 
ization. However, it exhibits blazar-like optical behaviour at 
low levels, diluted by thermal radiation (Impey et al. 1989, 
Valtaoja et al. 1990, Valtaoja et al. 1991) and has therefore 
been dubbed a mini-blazar by Impey et al. (1989). As expected 
within the context of the unification scheme, 3C 273 has a 
rather large viewing angle, estimated to be around 10 degrees 
(Lahteenmaki & Valtaoja 1999, Savolainen et al. 2006). This is 
smaller than for ordinary quasars, but larger than for blazars; 
for 3C 279 the estimated viewing angle is around 2 degrees 
(Lindfors et al. 2005 and references therein). 

In the simplest unification models, orientation is the only 
parameter. One might therefore expect that the borderline 
blazar 3C 273 would appear similar to 3C 279 if the view- 
ing angle were decreased. However, some of the differences 
we find between 3C 273 and 3C 279 (non-adiabatic jet flow 
in 3C 279, adiabatic in 3C 273; magnetic energy density lower 
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Fig. 2. Logarithmic evolution of the 
synchrotron emission of an average 
outburst in 3C 279. The three dimen- 
sional (log v, log S , log t) representa- 
tion is shown in (a), while the Cartesian 
projections of this surface are shown 
in (b), (c) and (d). b) Frequency ver- 
sus time representation with contours 
starting at a flux density of 0.02 Jy and 
spaced by 0.3 dex. c) The synchrotron 
spectra at different times spaced by 
0.2 dex. d) The lightcurves at differ- 
ent frequencies spaced by 0.2 dex. The 
thick solid line with arrows traces the 
time evolution of the spectral turnover, 
whereas the dashed line connects the 
peak of the lightcurves at different fre- 
quencies. The filled circles refer to the 
specific characteristics of the outbursts 
(see Section 4.2) and the vertical bars 
in (b) and (c) show the frequency cov- 



Frequency (GHz) erage of the 19 lightcurves. 




Time [ year ] 

Fig. 3. Model lightcurves of the average outburst in 3C 279 
at different frequencies spaced by 0.1 dex (grey lines). The 
six highlighted lightcurves are at frequencies, v, defined by 
log(v/GHz)= 3.5, 2. 5,. ..,1.0, in order of increasing timescales. 
t r and t p refer to transitions from Compton to synchrotron 
stage (t r =0.45 years) and from synchrotron to adiabatic stage 
(t p =3.34 years) respectively, tf refers to the time when the flat- 
tening of the spectral index starts (tf =0.046 years). 



than electron energy density in 3C 279 whereas in 3C 273 there 
is little difference) seem intrinsic to the sources and cannot 
be explained by differences in viewing angles alone. However, 



we note that variations in the viewing angle as a function of 
time/radial distance could play a role in the observed differ- 
ences between these two sources. 

4.2. Characteristics of Different Outbursts 

The evolution of each of the 15 individual outbursts is allowed 
to deviate in scale from the average outburst. For this we need 
three parameters corresponding to the logarithmic shifts along 
the three axes of Fig. 2, namely time, A log f, frequency, A log v 
and flux density, A log S . The values obtained for these shifts 
are given in Table 1 and their distribution is shown graphically 
with the filled circles in Fig. 2. 

We attempted to understand the physical origin of the dif- 
ferences seen between the outbursts by looking for correla- 
tions between A log S , A log v and A log t . Contrary to results 
for 3C 273 (Tiirler et al. 1999) we find no correlation between 
any of the shifts, and Fig. 2 shows clearly that the shifts do not 
align with any axis. It is also clear that the outbursts do not dif- 
fer mainly in amplitude as found for GRS 1915+105 (Tiirler et 
al. 2004) but also in peaking frequency and duration. 

To establish that our decomposition is not purely mathe- 
matical but corresponds to some physical reality, we compare 
the outbursts suggested by our fit to observed VLBI compo- 
nents. 3C 279 has been monitored with VLBI since the 1980s, 
and for the period considered in this paper the VLBI data cov- 
erage is good. In Table 2, the start times (To) of the outbursts, 
as derived from our modelling, are compared with the extrapo- 
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Shock nr. 


T 


A logS 


A logy 


A log? 


1 


1989.72 


0.25 


-0.18 


0.54 


2 


1990.57 


0.35 


0.03 


-0.45 


3 


1991.25 


-0.35 


0.48 


0.04 


4 


1992.51 


-0.28 


-0.57 


-0.03 


5 


1993.14 


0.18 


-0.25 


-0.29 


6 


1993.57 


0.10 


-0.08 


0.02 


7 


1993.98 


-0.06 


0.49 


0.03 


8 


1995.05 


-0.06 


-0.06 


-0.06 


9 


1995.86 


0.37 


-0.29 


-0.13 


10 


1996.09 


-0.38 


0.33 


0.20 


11 


1997.12 


0.37 


-0.70 


0.26 


12 


1997.14 


-0.12 


0.27 


0.22 


13 


1998.09 


-0.10 


0.23 


-0.24 


14 


1998.48 


0.04 


-0.21 


-0.10 


15 


1998.85 


-0.31 


0.41 


-0.07 



Table 1. The start times of the outbursts (r ) and their loga- 
rithmic shifts from the average shock evolution in flux density 
(A log S), frequency (A log v, negative if the outburst peaks at 
lower frequencies than the average outburst) and time scale (A 
log t, negative if the outburst evolution is faster than the evolu- 
tion of the average outburst). 



lated zero-epochs of the VLBI components (Wehrle et al. 2001, 
hereafter W0 1; Jorstad et al. 2004, hereafter J04). We also com- 
pare the observed flux densities of these knots (from W01 and 
Jorstad et al. 2005, hereafter J05) to the model flux densities 
of the outbursts (at the same date and frequency as the VLBI 
observations). We adopt the naming scheme of the VLBI knots 
from W01 (components C4 to C7a) and J04 (components C8 to 
C13). We note that there is a difference in identification of com- 
ponent C9 between these two papers and we choose to adopt 
the identification from the latter. 

Throughout the 1990s the VLBI maps showed two perma- 
nent features: C4 and C5. C4 is a very long-lived feature in 
3C 279. The flux density of this component has not decreased 
steadily with time but has dimmed or brightened on several oc- 
casions. This component contributes significantly to the initial 
flux density decay (which is a sum of all outbursts peaking be- 
fore 1989.0) but our model cannot account for its brightenings. 
For C5 there is no zero-separation epoch available because its 
fitted motion is consistent with zero speed. The flux density of 
this component was falling throughout the 1990s and it disap- 
peared in 1998. We suggest that the first, very long-lived out- 
burst in our fit is related to this stationary component and is not 
caused by a moving knot. 

We identify the other components by looking at the closest 
match between the zero-separation epoch of the VLBI compo- 
nent and the start time of the outburst. We find that for most 
of the components the match is rather good, although in most 
cases To is outside the error-bars of the extrapolated VLBI 
zero-separation time. There are large mismatches between the 
observed VLBI component flux densities (from W01 and J05) 
and the flux densities derived for our model components. The 
model outbursts have systematically higher flux densities than 
the real VLBI components. On the other hand, the flux density 
of the underlying jet component in our model is lower than the 



flux density of the VLBI core at any given epoch. This sug- 
gests that the underlying jet component is underestimated in 
our decomposition, while the flux density of the outbursts is 
overestimated during their decay stage. 

We can see from Table 2 that there are four outbursts 
without a corresponding VLBI component. The third and the 
seventh outbursts are high-frequency peaking outbursts (see 
Table 1) which are probably too weak to be seen in the maps 
by the time the decaying shocks are emitting at VLBI frequen- 
cies. On the other hand, the sixth outburst should be seen in the 
VLBI map but when we look at the 22 and 37 GHz lightcurves 
(those closest to the VLBI frequencies) we can identify the 
fifth and the sixth outburst as a single event. The millimeter 
lightcurves cannot be fitted with one outburst and we therefore 
suggest that in higher frequency VLBI we might see two com- 
ponents. Another possible explanation is that the fifth outburst 
is just a core flare (it has rather short duration, see Table 1) and 
that the VLBI component should be recognized with the sixth 
outburst. The VLBI observations of W01 show that the core 
flux density increased during 1993, thus supporting the latter 
scenario. The long-lived tenth outburst, which is rather faint 
at low frequencies, is likely to be related to the ninth outburst 
and to the exceptionally bright VLBI component C8. In fact, it 
might be the case that we see the tenth outburst of the fit be- 
cause of the brightening of the component C8 seen in VLBI 
maps by W01 and J04. The mismatch between the flux densi- 
ties of the outbursts and VLBI components is very severe dur- 
ing 1997 (eleventh, twelfth and thirteenth outbursts). VLBI ob- 
servations from this period suggest that the direction of the jet 
"nozzle" changed (J04). Our simple model cannot reproduce 
this behaviour and therefore the discrepancy in component flux 
densities is bigger than for other outbursts. 

The fifteenth outburst peaks at a high frequency, similar to 
the third and the seventh outburst, therefore it has a very low 
flux density at VLBI frequencies and should not be visible in 
the maps. It is therefore likely to be unrelated to the VLBI com- 
ponent CI 3, which then has no outburst counterpart in our fit. 
This is probably because the lightcurves we use in our fit end 
soon after the zero epoch of CI 3. As the last flux density points 
at 22 and 37 GHz are decreasing we cannot fit an additional 
outburst there even if one was beginning at that time. 

To make a further comparison of our model fit to VLBI ob- 
servations we produced artificial time-series of VLBI maps of 
3C 279 at 22 GHz. These five maps describing the evolution 
of our model jet are shown in Fig. 4. The flux density of the 
underlying jet component of our fit is fixed to the core and re- 
mains constant at a level of 5.8 Jy. The components of the map 
correspond to model outbursts of our fit and we assume them 
to move with a constant apparent velocity of 5c (W01 found 
that components have apparent velocities of 4.8c to 7.5c). The 
maps are scaled to the distance of 3C 279 assuming Ho = 70 
km s _1 Mpc -1 , and their angular resolution is set to 0.2 mas as 
in Fig. 3 of W01 to ease the comparison. 

At first sight, the simulated maps look very different to the 
observed ones. This is, however, mainly due to the fact that the 
components starting before 1989 could not be included in the 
simulated maps and explains the growth of the simulated jet as 
well as the absence of component C4. Another clear difference 
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2.86* 


1996.02 
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9 
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C8 
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1996.02 
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1996.09 
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1997.12 


C9 
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2.959 


1998.23 


43 


0.81 


12 
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C10 


1997.42+0.1 


1.0 
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1998.23 


43 


3.21 


13 
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Cll 


1997.59+0.1 


0.4 


4.305 
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43 


1.02 


14 


1998.48 


C12 


1998.50+0.1 


5.2 
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1998.94 


43 


0.51 


15 


1998.85 


C13 


1998.98+0.1 


0.05 


2.295 
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43 


2.16 



*sum of C7a and C8 

Table 2. Comparison of the outbursts found here with observed VLBI components. The start times, Tq, of the outbursts are 
compared to the zero-separation epochs, 7o(knot), of VLBI components derived by Wehrle et al. (2001) and Jorstad et al. (2004). 
We also compare the extrapolated flux density, F exp , of the outbursts at the time of the VLBI maps to the observed VLBI flux 
density, F D b s (Wehrle et al. 2001; Jorstad et al. 2005). The ninth column gives the rise times of the outbursts at 43 GHz. 



is the bright outer blob in the simulated maps, which corre- 
sponds to the first, very long-lived outburst. As stated above, 
this is likely not a real propagating shock wave, but accounts 
for the flux density of the stationary component C5 and, maybe 
also to some extent, to the strong long-lived and re-brightening 
component C4. Apart from this, we note a predominance of 
the core emission in the observed VLBI maps compared to our 
simulated maps. This illustrates the fact that the flux density of 
the core in the VLBI maps of W01 is higher at all times than 
the flux density of our underlying jet component, and it varies 
from 8.8 to 21.1 Jy. W01 found that flux density variations of 
the core are largely responsible for the total flux density vari- 
ability observed by single-dish antennae. The core brightenings 
are often followed by the ejection of superluminal components. 
The model maps presented here suggest that during the initial 
stage of the model outbursts it is not possible to resolve them 
from the VLBI core; they are first seen as a brightening of the 
core. In the later stages, the longer lasting outbursts are seen 
as VLBI components while short-lived components might only 
be seen as core-flares due to the limited resolution of the instru- 
mentation. 

Using the shifts given in Table 1, we can calculate the rise 
time of each outburst at 43 GHz, defined as the duration from 
the beginning of the flare to the time when the maximum flux 
density at 43 GHz is reached. From the values given in Table 2, 
we note that most of the outbursts peak within 1 year and they 
have a maximum flux density amplitude at rather low frequen- 
cies. Using the same assumptions as above, the proper mo- 
tion is 0.16 mas/yr, so most of the flares peak while they are 
still within the core (even if we would assume a resolution of 
> 0.1 mas like in J05). When the component is resolved in 
VLBI maps, the flux density is already decaying but the time 



elapsed since the shock's maximum stage is short and the com- 
ponent is still well detectable with VLBI. On the other hand, 
the outbursts with a very long rise time at 43 GHz are associ- 
ated with shocks peaking at high frequencies which are heavily 
self-absorbed at VLBI frequencies. After several years, when 
the self-absorption turnover reaches 43 GHz, the component 
flux density has already decreased appreciably and often re- 
mains below the VLBI detection limits. 

We conclude that the outbursts of our fit correspond sat- 
isfactorily to real observed components in the jet. This is in 
accordance with previous studies: Tiirler et al. (1999) found 
that their model outbursts correspond to VLBI components 
in 3C 273, and Savolainen et al. (2002) studied a sample of 
27 AGN and found a clear connection between the observed 
VLBI components and the total flux density outbursts in 22 
and 37 GHz radio lightcurves. 

4.3. Comparison with EGRET gamma-ray states 

As 3C 279 is the most frequently observed EGRET blazar it is 
very interesting to take a closer look at the components present 
at different EGRET epochs. It is reasonable to assume that most 
of the gamma-ray emission will originate in the newest shock 
component since, during the early stages of shock evolution, 
the dominant cooling mechanism is inverse Compton losses. 
This shock is usually also the dominant component in the 
submm to optical range of the spectrum (see below). Therefore 
we calculated the distances of these components along the jet 
with: 

L * D 2 ■ c{\ + zT x M obs (1) 

where z=0.538 for 3C 279 and for simplicity we assume D = 
10 for all shock components, in accordance with Lindfors et al. 
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Fig. 4. A time sequence of five model VLBI maps of 3C 279 
at 22 GHz between 1990 and 2000. As in Fig. 3 of W01, the 
map resolution is 0.2 mas, with a lowest contour at 50 mJy 
beam -1 and subsequent contours increasing by factors of 2. 
Each component (marked with x) is assumed to move with a 
constant apparent velocity of 5c with respect to the core of the 
jet (marked with +). An animation of this figure is shown at: 
http://isdc.unige.ch/~turler/jets/ 
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A?,,;,., 
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L[pc] 


P5b (1996.092) 


10 


0.006 


very large flare 


0.12 


P8 (1999.070) 


15 


0.195 


high 


3.88 


P5a (1996.063) 


9 


0.206 


high 


4.10 


PI (1991.47) 


3 


0.225 


high 


4.48 


P3a (1993.858) 


6 


0.288 


moderate 


5.74 


P6b (1997.470) 


12 


0.328 


moderate 


6.53 


P3b (1993.979) 


6 


0.409 


moderate 


8.15 


P2 (1993.004) 


4 


0.49 


low 


9.76 


P6a (1997.010) 


10 


0.924 


low 


18.41 


P4 (1994.970) 


7 


0.99 


low 


19.72 



Table 3. The EGRET epochs ordered by increasing A? t, s 
[years], the time interval since the start of the last outburst as 
derived here. This order is compared to the gamma-ray state 
adopted from Hartman et al. (2001). Also an estimate of the 
distance, L, of these shock components from the apex of the jet 
is given (see text). 



(2005). Af b s is the difference between the onset time of the out- 
burst (To) and the time of the EGRET observation. The results 
are shown in Table 3. 

We note that all but one (outburst 10 at epoch 5b) of the 
shock components are far beyond the broad-line region (BLR) 



boundary of 0.4 pc (Kaspi et al.2000) at the time of the EGRET 
observation. Therefore, if the gamma-ray emission is indeed 
related to shocks travelling in the jet then mechanisms based on 
external photon fields reprocessed in the BLR - actually most 
of the External Compton mechanisms suggested up to date (e.g. 
Sikora et al. 1994; Hartman et al. 2001, and references therein) 
- are ruled out for all other epochs but 5b, when a very large 
flare in gamma-rays was observed. 

There is a clear pattern in Table 3 when we look at the 
Af b s and compare them to the gamma-ray state. The high states 
(and the very large flare in epoch 5b) correspond to the shortest 
Afobs, while the low states correspond to the longest Af b s and 
the moderate states fall neatly between these two cases. This 
pattern clearly suggests that the gamma-ray emission is indeed 
connected to the shocks travelling in the jet and that the high- 
est gamma-ray fluxes are observed when there is an early stage 
shock present in the jet. 

To gain confidence in the reality of this pattern, we also 
looked for a correlation with the next outburst beginning after, 
rather than before, a gamma-ray epoch. But in this case, there 
is apparently no pattern between the gamma-ray state and the 
time delay between the gamma-ray epoch and the beginning of 
the outburst. 

It is instructive to compute the probability that the 10 Af b s 
are grouped as they appear in Table 3 by chance, assuming no 
real correlation exists between Af b s and gamma-ray state. We 
simulate 100,000 sequences of the 10 Af b s and repeated the 
simulation 1000 times. We found that in average the Af D b s are 
ordered as in Table 3 24.45 times per simulation, which cor- 
responds to a probability of 0.024% or 3.7<x away from the 
excepted value (assuming a normal distribution). 

In the simulations we assume that Af Q b s and the gamma- 
ray state are well defined quantities with no uncertainties. For 
low -peaking outbursts, the error in the determination of To (and 
therefore Af b s ) is rather small because the radio lightcurves are 
very well sampled. For high-peaking outbursts the uncertainty 
is larger so, for example, we cannot exclude the possibility that 
outburst 10 would start 0.01 years later and the dominating out- 
burst at epoch 5b would be outburst 9 with Af D b s = 0.235. This 
would still not destroy the pattern seen in Table 3, but demon- 
strates that although the observed correlation has a high statis- 
tical significance, it should still be considered suggestive rather 
than definite. 

It has already been established statistically for a sample 
of sources that gamma-ray flares are connected with growing 
shock components propagating in the jet (Jorstad et al. 2001, 
Lahteenmaki & Valtaoja 2003). Our finding supports this sce- 
nario and suggests additionally that there is a dependence be- 
tween the gamma-ray flux and the distance of the shock compo- 
nent from the apex of the jet. We emphasize the need of further 
studies on this issue when future experiments provide gamma- 
ray lightcurves with better time coverage. 

To discuss further the link between gamma-ray emission 
and synchrotron outbursts we show in Fig. 5 the low-frequency 
spectral decomposition at each EGRET epoch as derived from 
our lightcurve fitting. For each epoch, the total model spectrum 
and its decomposition into individual outbursts is compared to 
the average spectrum of 3C 279 (points). We see that at two 
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Fig. 5. The synchrotron spectra of 
3C 279 at 10 EGRET epochs. The 
points indicate the observed mean flux 
density levels. The solid line is the total 
model spectrum at each epoch which is 
the sum of the spectra of the underlying 
jet component (dashed line), the overall 
spectrum of the outbursts peaking be- 
fore 1989.0 (thick dotted line) and the 
individual outbursts (dotted lines). 



of the high gamma-ray state epochs (PI and P8) the model 
spectrum in the optical range is clearly above the average ob- 
served flux density levels, but at two other epochs (5a and 5b) 
the model spectrum fits very well the average optical flux den- 
sity. At these epochs it is instead the mm/submm flux density 
which is above the average, which is not seen in epochs PI and 
P8. This is easy to understand by looking at the dominant shock 
components at these epochs: at PI and P8 the last outburst is 
high-frequency peaking, while the ninth outburst at epochs 5a 
and 5b peaks in the millimeter range. 

There are four epochs when the fitted optical spectrum is 
below the observed mean flux density level. Three of these 
epochs correspond to a low gamma-ray state and the fourth to 
a moderate gamma-ray flux (epoch P6b). There are also two 
epochs (P2 and P6a) corresponding to low gamma-ray states 
and for which the millimeter to submillimeter spectrum is be- 
low the observed mean spectral energy distribution. Finally, we 
note that the centimeter spectrum is dominated by the underly- 
ing jet component at all epochs, and that the fitted spectrum is 
below the observed one at the first two epochs but above it at 
the last three epochs. 

In summary, the pattern we see in Fig. 5 is in accordance 
with and strengthens the hypothesis that gamma-ray emission 
originates in young, growing shocks which are dominant in the 
mm-to-optical regime at the time of high gamma-ray states. 



5. Conclusions 

We have studied synchrotron flaring in the jet of 3C 279 by 
decomposing the multifrequency lightcurves into a series of 
outbursts. The analysis is performed within the generally ac- 
cepted shock-in-jet scenario, with the evolution of the model 
outbursts following the Marscher & Gear (1985) prescrip- 
tion. Similar methods have been used previously to study the 
quasar 3C 273 (Tiirler et al. 1999, 2000) and the micro-quasar 
GRS 1915+105. For 3C 279 we find: 

1 . The lightcurve decomposition supports the Marscher & 
Gear model. The smooth low-frequency lightcurves are espe- 
cially well reproduced by the model. In the high-frequency 
(millimeter to optical) domain, the sharp peaks of the observed 
lightcurves are less well reproduced by the model lightcurves 
during the strongest outbursts. Although the model used here 
makes a lot of assumptions and simplifications which might 
not hold in detail, it is still clear that there is some mecha- 
nism at work which is not well described by the Marscher & 
Gear (1985) model, and which produces the sharp peaks of 
outbursts at high frequencies. This is in accordance with previ- 
ous findings (e.g. Valtaoja et al. 1999): both the theoretical and 
numerical jet models always produce smooth peaks in model 
lightcurves while the observed outbursts have exponential rises, 
sharp peaks and exponential decays. 

2. The outburst evolution in 3C 279 follows the character- 
istic three stage pattern but instead of the synchrotron plateau 
stage seen in 3C 273, the flux density is already decreasing dur- 
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ing the second stage. This might be explained by a highly non- 
adiabatic jet flow, where radiative losses play a significant role. 
We also find that the lightcurve of an average outburst is dif- 
ferent for 3C 273 and 3C 279. The differences found between 
these two sources appear to be mainly intrinsic and not due to 
different orientations. The borderline blazar 3C 273 is thus not 
simply 3C 279 seen from a larger viewing angle. 

3. We find a rather good agreement between the zero- 
separation epoch of VLBI components and the start times of 
model outbursts. However, some high-frequency peaking out- 
bursts do not have a corresponding VLBI component which can 
be explained by the long delay before their spectral turnover 
reaches VLBI frequencies. 

4. When comparing the gamma-ray state at different epochs 
from Hartman et al. (2001) to the time elapsed since the start 
of the last synchrotron outburst we find that the shortest time 
intervals (corresponding to the shortest distances from the apex 
of the jet to the shock) are seen during the highest gamma-ray 
states, and the longest ones during the lowest gamma-ray states. 
This supports previous findings (e.g. Valtaoja & Terasranta 
1996; Jorstad et al. 2001; Lahteenmaki & Valtaoja 2003) that 
the gamma-ray flares are connected to shocks propagating in 
the jet. Although during the highest states the shock compo- 
nents are rather new, they are still well outside the outer limits 
of the broad-line region. This suggests that external Compton 
models based on the soft-photon field originating from the BLR 
are unlikely to be able to explain the gamma-ray flares in this 
source. 
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